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. Abstract 

^ I It has been shown recently that Bose Gase with weak pair (enough 

■ weh) interaction is non ergodic system. But Bose Gase with weak pair 

^ , interaction is so general system that it is evident that the majority 

Q I of statistical mechanics systems are non ergodic too. It is also has 

been shown that it is possible to generalize the scheme of standard 
statistical mechanics and thermodynamics to take into account non 
' ergodicity. This generalization is called a generalized thermodynam- 

' ics. In some points this generalized thermodynamics coincide with 

. standard equilibrium thermodynamics but some new specific results 

^ ; take pi ace. It has been shown that this new generalized thermody- 

namics can be used to explain some physiological phenomena which 
take place in the living cell when the cell is exciting and dying. 

In the present paper we try to illustrate some basic points of this 
generalized thermodynamics on some physical examples. 
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1 Introduction 



The purpose of this paper to give some examples on so called generalized 
thermodynamics recently developed in [1,2]. 

In [1] it has been shown that even most reahstic systems of statistical 
mechanics (for example usual Bose Gase with weak (enough well) pair inter- 
action) are non ergodic systems. These means that there exist non trivial 
first integrals of the system commuting with momenta and number particle 
operatots. It has been recently shown [2] that it is possible to generalize the 
scheme of standard statistical mechanics and thermodynamics to take into 
account the existence of such non trivial first integrals. These generalization 
we call the generalized thermodynamics. It has been also shown that this 
generalized thermodynamics is very useful to explanation of some physiolog- 
ical phenomena which take place in the living cell when the cell is excited 
dying. 

Let us describe the basic elements of the scheme of this generalized ther- 
modynamics. The non ergodic theorem states that for a wide class of re- 
alistic systems of statistical mechanics there exists non trivial commuting 
(in involution) first integrals Ki,...,Kn, N = 1,2,... commuting with mo- 
menta and number particle operators. In purpose of simplicity we will talk 
about Hamiltonian instead of Hamiltonian, momenta and number particle 
operators. The starting point of generalized thermodynamics is a following 
expression for distribution function (density matrix in quantum case): 



where a; is a point of phase space of the system and E, K[, ...^K'^ are the 
observable values of energy and integrals Ki^...,Kn- The entropy corre- 
sponding to this distribution is defined as a logarithm of statistical weight. 
The statistical weigh, by definition, is a number of microscopic configuration 
of the system, corresponding to a given macroscopic state i.e.: 



W{E,K[,...,K'^)= / dTJ{H{x)-E)Y[5{Ki{x)-K'^) (2) 



N 



p{x) = const d{H{x) - E)Y[6{Ki{x) - K-), 



(1) 



N 



i=l 



and 



(3) 
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where dT^ is an element of phase volume. The expression ([T]), and 
are the generalization of standard microcanonical Gibbs distribution and its 
entropy to the case when there exists non trivial commuting first integrals 
Ki, Km of the system. 

Let us explain why we require pairwise commutativity of integrals 
Ki,...,Kiy or why we require that Wi,j = 1,...,N {Ki,Kj) = where (■,■) 
denotes the commutator in quantum case or the Poison bracket in classical 
case. We want to describe the state with definite values K[, K'^ of first 
integrals i^i, Kj^. So the integrals Ki, i = 1, N must be contemporane- 
ously measurable. But in quantum mechanics this means that the integrals 
Ki,...,Ki\f must be pairwise commutative. The requirement that the inte- 
grals Ki, ...,K]\r must be in involution (in classical case) is clear now from 
the remark that the Poison bracket is a classical analog of commutator. 

Now let answer the question when generalized thermodynamics gives re- 
sults which differs from results following from standard thermodynamics. 
I.e. let us answer the question when using of distribution ([1]) leads to results 
which differ from result obtained by using standard microcanonical Gibbs 
distribution. 

For simplicity consider the case = 1. The case of an arbitrary can 
be considered by analogy to this case. For a fixed E S{E, K') is a function 
of K' . There exist two typical cases of behavior of this function. 

1) S{E, K') (for fixed E) has a maximum in isolated point K' = K". 

2) S{E, K') achieve a maximum at whole interval K' G [a, h] of nonzero 
length. 

There arise a question, which observable values K' of integral K could be 
realized in nature. It is essentially to suppose that it could be realized such 
and only such values K' of K which corresponds to the maximum of entropy 
(for fixed energy). This is one of main point of generalized thermodynamics. 
It has been shown [2] that in case 2) generalized thermodynamics gives no 
more than standard equilibrium thermodynamics {K') is uniquely defined). 
But new interesting physics arise in case 2). Note that the situation when 
thermodynamical functions have a plateau as functions of their arguments is 
a typical situation in the theory of phase transitions. 

Some general results of standard equilibrium thermodynamics take place 
in generalized thermodynamics [2] . For example a well-known relation 

dE = TdS - PdV (4) 

takes place in new situation. Here E is an energy of the system, T is a 
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temperature of the system, K is a volume of the system, P is a pressure. 
But some new specific resuhs take place for generalized thermodynamics. 
These new results have been used [2,3] for explanation of several physiological 
phenomena which takes place in biological cell when the cell is dying or 
moving from resting state into activated state. 

In the present paper we give some new physical examples to illustrate 
some basic points of this generalized thermodynamics. The paper is com- 
posed as follows. In section 3 we describe some new derivation of the Gibbs 
distribution from the property of asymptotical factorization of correlations 
based on ideas of nonequilibrium renormalization theory. In section 3 we 
reformulate the Bardeen — Cooper — Schrieffer model on the language of 
the generalized thermodynamics. In section 5 we describe some derivation of 
the Boltzmann kinetic equation based on some ideas of N.N. Bogoliubov 
concerned with thermalization in oscillator interacting with thermostat. In 
section 6 we describe how to obtain non-trivial first integrals for the equa- 
tions on order parameter of the system in many-phase domain in the theory 
of second order phase transitions. Section 6 is a conclusion. 

2 Derivation of Gibbs distribution from the 
property of asymptotic space factorization 
of correlations. 

In the present section we prove that all stationary translation-invariant states 
of Bose Gase with weak (enough good) pair interaction satisfying to the 
property of (enough fast) asymptotic space factorization of correlations are 
the Gibbs states (described by the Gibbs distribution). 

But at first let us recall the main result of the paper [1] and demonstrate, 
how the non ergodic property follows from this result. 

The main result of [1] can be formulated as follows: 

Theorem 1. For Bose Gase with weak pair interaction with kernels 
from Schwartz space in the sense of formal power series on coupling constant 
there exists non- Gibbs functional (•), commuting with the number of particle 
operator such that the correlators 

are translation invariant, do not depend on t and satisfy the weak cluster 
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property. Here ^ are the secondary quantized wave function and its complex 
conjugated ^"(x) := ^(x) and the weak cluster property means the following 



lim / {^^{t,xi + 6ieia)...^^{t,Xn + 5neia))f{xi, 

|a|->oo J 

= J (*±(t,XiJ...*±(t,a;,J)(*±(i,a:,J...*±(i,a:,J) x/(xi, 
there 5i e {1, 0}, i — 1, 2...n and 

ii <i2 < ... < ifc, 

ik+l < 'i-k+2 < ■■■ < in, 

{ii,i2, ik} = {i = l, 2...n\5i = 0} 7^ 0, 
{ik+i, ik+2, in} = {i = l, 2...n\6i = 1} 7^ 0. 

f{xi, Xn) is a test function (i.e. the function from the Schwartz space), ei 
is a unit vector parallel to the x-axis. About secondary quantization see for 

example [4, 5]. 

Let us prove that the existence of such functionals implies non-ergodic 
property of the system. The more accurate proof of this fact see in section 
10 of [1]. Suppose that our system is ergodic, i.e. there are no first integrals 
of the system except energy. Then, the density matrix p of the system 
corresponding to the functional (■) is a function of energy. We can represent 
this density matrix p as follows: 

where H is a. Hamiltonian of the system and the sum can be continuous 
(integral) . Let 1 be some enough large but finite subsystem of our system. Let 
2 be a subsystem obtained from 1 by translation on the vector / of sufficiently 
large length parallel to the Let 12 be a union of the subsystems 1 

and 2. Let pi, p2 and pi2 be the density matrix of the subsystems 1, 2 and 
12 respectively. By the same method as the method used for the derivation 
of the Gibbs distribution we find: 

E e e 



••;^?7.)^ X\...d Xfi 
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in the obvious notation. Here Hi and H2 be the Hamiltonians of subsystems 
1 and 2 respectively. But the weak cluster property implies that 

P12 = Pi ® P2- 

Therefore all the coefficients are equal to zero except one. We find that 

p = c5{H - Eo) 

for some constants c and Eq. So each finite subsystem of our system can be 
described by Gibbs formula and we obtain a contradiction. 

But the state (not necessary positive defined) mentioned in theorem 1 
satisfies only to the weak cluster property i.e. to the property of asymptotic 
space factorization of correlations only " in one direction" . 

It follows from the proof of theorem 1 from [1] that we can achieve, that 
(•) will satisfies to the property of asymptotic factorization of correlations 
"in two directions". But there arise principled difficulties if one try to prove 
that the state (■) satisfy to cluster property in "all (three) directions". 

It has been mentioned above that all stationary translation-invariant 
states satisfying to (enough fast) cluster property are the Gibbs states. Let 
us now recall the standard derivation of this fact. Then we will point out 
some problems connected with this derivation and present our new proof of 
this fact based on results of our renormalization theory of nonequilibrium 
(Keldysh) diagram technique. 

Let us at last describe the standard derivation of Gibbs distribution (from 
cluster property) at classical level [6]. This derivation is based on the fact 
that there no exist additive first integrals of the system (linear) independent 
of Hamiltonian H, momenta P and angular momenta M. Let p be a station- 
ary translaion-invariant distribution function of the system satisfying to the 
cluster property of the system. If the system is divided into two subsystems 
1 and 2 cluster property implies that 

P = P1P2, (5) 

where pi and p2 are distribution functions for subsystems 1 and 2 respectively. 
In other word the distributions for subsystems 1 and 2 are independent. 
Therefore In p is an additive integral of motion and can be represented as 
linear function of H, P, M. In other words 

In p^aH + ^P + ^M (6) 
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for some real number a and real vectors /3 and 7. Therefore 

But we assumed that p is translation-invariant. Therefore 7 = and 

p = e"^+^^^ (8) 

The distribution ([8]) is a standard Gibbs distribution. But there are two 
problems with this derivation. The first one, this derivation is performed 
only at classical level. The second one, we could not find anywhere the proof 
of the fact that Hamiltonian H momenta P and angular momenta M are the 
complete linear independent set of additive first integrals. 

Let us present now our derivation of Gibbs distribution (from the cluster 
property). This derivation is not rigorous too but we hope it is of some 
interest. Note that this derivation uses some basic ideas of Haag-Ruelle 
scattering theory [7]. 

But at first let us give some previous definitions. 

Definition 1 .Let S'(M^) be a Schwatrz space of test functions (infinitely- 
differentiable functions decaying at infinity faster than any inverse polyno- 
mial with all its derivatives). The algebra of canonical commutative relations 
C is an unital algebra generated by symbols a~^{f) and a(/) / e S'(M'^) sat- 
isfying the following canonical commutative relations: 

a) a^if) is a linear functional of /, 

b) a{f) is an antilinear functional of /, 



[a(/),a(<7)] = [a+(/),a+(/)] = 0, 
Hf),a+{g)] = {f,g), 

where {f,g) is a standard scalar product in L^(]R^), 

{f,g) := J r{x)g{x)d'x 

Remark. We will widely use generalized "elements" of the algebra of canon- 
ical commutative relations a{k), a'^{k) defined according to the following 
relations 



a{f) = J a{k)f{k)d% 
a^if) = [ a+{k)t{k)d^k. (9) 
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a{k), a'^{k') are called the annihilation-creation operators and satisfy to 
the following canonical commutative relations: 



[a+{k),a+{k')] = [a{k),a{k')] = 0, 
[a{k),a{k')] = 5{k-k'). 



(10) 



Definition 2. The field operators ^(x), ^"•"(x) are defined as follows 



The rigorous definition of ^(x) and ^"'"(a;) could be obtained from this defi- 
nition by using the notion of the Fourier transform of distributions. 

Definition 3. Let p be a state on the algebra of canonical commuta- 
tive relations (the algebra generated by smoothed secondary quantized wave 
functions) (CCR-algebra). We say that p is a Gauss state if we can calcu- 
late its values at elements of CCR-algebra by using the Wick (Bloch — De 
Dominicis) theorem through pair correlations. 

Now let us recall some notions connected with nonequilibrium (Keldysh) 
diagram technique. Let p be an arbitrary state on the algebra of canonical 
commutative relations. Let us introduce the Green functions for the system 
(corresponding to this state). 



Here ti,...,tn are times, symbol H near "if^ means here that are 
operators in Heizenberg representation and the symbol T is a symbol of 
chronological ordering. 

We will consider the system describing by the following Hamiltonian 




(11) 



p{T{^%{ti,Xi),...,^%{tn,x„))). 



H^Ho + XV, 



(12) 



where Hq is a free Hamiltonian 




(13) 
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oj{k) = ^. (14) 
A e R is a coupling constant and 

y = ^ y d^x(fx'<i/+{x)<i/+{x')V{x - x')*(x')*(x), (15) 



where V{x) is an arbitrary test function from S'(M ). In nonequilibrium di- 
agram technique we require the following representation for the Green func- 
tions: 

p(T(^|(ti,xi),...,^|(t„,x„))) = 

Po(^-'r(*±(tl, Xi), *±(t„, Xr^)S)). (16) 

The symbol near means here that are operators in the Dirac rep- 
resentation (representation of interaction). The 5- matrix has the form 



H-oo 

5 = rexp(-i J V^{t)dt), 



and 

+ 00 



S-^ = fexp{i j V„{t)dt). 



Here T and T are symbols of the chronological and the antichronological 
ordering respectively, po is some Gauss state defined by density function 
n{k) as follows 

Po{a+{k')a+{k)) = po{a{k)a{k')) = 0, 

Po(a+(k)a(k')) = n(k)5(k - k'). (17) 

The state po is called the asymptotical state. 

Let us recall the basic elements of nonequilibrium diagram technique. 
The vertices coming from T-exponent are marked by symbol — . The vertices 
coming from T-exponent are marked by symbol -|-. There exist four types of 
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propagators 

G^'iti - t2, xi - X2) := po(*(ii, a^i)*+(i2, ^2)), 

GQ-^ih - t2, Xi - X2) := Po(*+(t2, X2)*(tl, Xi)), 

Go~{h - t2, xi - X2) := Po(T(*(ti, xi)*+(t2, 3^2))), 

^^+(^1 - t2, Xi - X2) := Po(T(*(tl, Xi)*+(t2, X2))). 

Let us write the table of propagators 

r--(f r)-i f ^'fe r l + ^(^) i.-^(c..-M 

° ^' ^" 7 (2^)4^0; -a;(A;) + zO a; - c^(A;) - zO^ 



ci^A; 1 + n(A;) 



{ .. . , - .. \.n K 



-i{wt—kx) 



{2tiY^u - u{k) + iQ uj-u{k)-iO- 

It has been shown in [1] that (usually) there exists divergences in Keldysh 
diagram technique. A typical example of divergent diagrams is pictured at 
fig.l 

fig. 1 



• • • 



+v y + +v y + 



The ovals represent the sum of one-particle irreducible diagrams. These 
diagrams are called chain diagrams. Let us suppose that all divergences of 
self-energy parts (ovals) are subtracted. The divergences arise from the fact 
that singular supports of propagators coincide. 

For all ti, ^2 e let us define the evolution operator ^(ti, ^2) as follows: 

S{ti,t2) = e^t^Ho^-i{t^-t2)H^-itHo_ (^8) 

Note that, for example, that 

o 

5(0, -00) = rexp(-i J Vo{t)dt), 

—00 
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-l-oo 

5(+oo,0) = Texp(-i j Vo{t)dt) (19) 

o 

One can prove that the following lemma holds. 

Lemma 1. Vt G M the following equalities holds: 

e-''^S{0, -oo) = S{0, -oo)e-^*^o 

(20) 

and 

S{+oo, 0)e-**^ = e-''^'S{+oo, 0). (21) 
Now let us consider the equilibrium Green functions defined as follows 

^tr{T{'Ir^{t,,x,)...^^{tn,Xn))e-'^), (22) 



where 



H-uN , 

Z = tr(e — —), (23) 



T > is a temperature, /u G M is a chemical potential and is a number 
particle operator 

N = j a^{k)a{k)d?k (24) 
In other words we put in ([T6l) 

P(-) = |tr((-)e-^). (25) 

It follows from the lemma 1 that the following lemma holds. 
Lemma 2. 

= Po(^"'n*o^(ti,a:i),...,vi/±(t„,x„)5)), (26) 

where po is a Gauss state on the CCR-algebra defined as usual by its density 
function n{k) of the form: 

<k) = (27) 

e T — 1 
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It follows from physical reasonings that the equilibrium Green functions does 
not contain divergences. So we have the following lemma: 

Lemma 3. If the asymptotical state po in Keldysh diagram technique 
corresponds to the density function n{k) of Bose-Einstein form 

n{k) = (28) 

e T — 1 

then the Keldysh diagram technique does not contain divergences and cor- 
responding green function are equilibrium 

p{T{^f%{h,Xi),...,^%{tn,Xn))) 

= itr(T(!P'±(t„ x,)...^^{tn, Xn))e-^), (29) 

Now let us prove the fact that if the Keldysh diagram technique does not 
contain divergences then the corresponding Green functions are equilibrium, 
i.e. the asymptotical state po has Bose-Einstein form. In order to do this let 
us analyze the divergences of Keldysh diagram technique in lowest possible 
order in A. Such divergences may come only from the diagrams pictured at 
fig. 2. 

fig. 2 




where ovals means here the sum 
of all self-energy diagrams of lowest possible order in A. It has been shown 
in [1] that these diagrams contain divergences for density function n{k) of 
general form. It has been also shown in [1] that these diagrams could be 
subtracted by the following renormalization of the asymptotical state 

1 - / ho{t)dt 

po(-) ^ ^Po(e — (•)), 

where 

h = [ h{k)a+{k)a{k)d?k, 
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h{k) is a real-valued function, 

ho{t) := e''^"he-'^"° (30) 

and 

- / h(t)dt 

Z^poie - ) 

for suitable h — J h{k)a^{k)a{k)(Pk. 

It also could be extracted from [1] that the divergences in considered 
diagrams exists if and only ii h ^ 0. But it has been also shown in [1] that 

- 2.(p)(l':nl)) ^'W- 

where St{p) is a collision integral (in lowest possible order of perturbation 
theory) and St{p) = if and only if n{p) has a Bose-Einstein form (see for 
example [8]). So we have proved the following 

Lemma 4. If Keldysh diagram technique does not contain divergences 
the the asimptotical state has a Bose-Einstein form 

Po{a+{k')a+{k)) = po{a{k)a{k')) = 0, 
Po{a+{k)a{k')) - n{k)6{k - k'), 

n{k) = (32) 

e T — 1 

and the corresponding Green function are equilibrium. 

Now let (•) be some stationary translation invariant state on the CCR- 
algebra, commuting with the number particle operator. Suppose that (•) 
satisfy to the asymptotic property of (enough fast) space factorization of 
correlations. Put by definition 

W = S(0,-oo). (33) 

Let (•)o be a state on the CCR-algebra such that 

(•) = {W-\-)W)o. (34) 

One has 

{■)o^{W{-)W-')- (35) 
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Let us formulate now our main lemma. 

Lemma 5. If the state (■) satisfies to the property of asymptotic (enough 
fast) factorization of correlations then (•)o is finite (does not contain diver- 
gences) and Gauss. 

Proof. We have 

lim e+^*^e-^*^° (36) 

t— >— oo 

and 

W-^^ lim e+'*^°e-'^". (37) 

These two equations could be simply established in the sense of formal series 
of perturbation theory and this enough for our aims. 
Therefore we have 

{■)o = {W{-)W-') 

— lim (^^+itiH^-itiHo^_^^+it2Ho^-it2H\^ 

ti— >~ oo, t2— >— oo 

= lim (e+^*^e-^*^°(-)e+^*^°e-^*^) 

t— >— oo 

= lim (e-'*^°(-)e+'*^°). (38) 

t—^—oo 

Last equality takes place because the state (•) is a stationary state. In result 

(•)o= lim (e-'*^<'(-)e+'*-^<'). (39) 

t— oo 

Let us now introduce so-called Wightman functions (distributions) 

■^Ij ^Ij •••) -^n; 

(7„), t eR xi, ...,Xn e B?, (Ji, ...,<7n G {+ , -} as foUows: 

W„,i(xi,(7i,...,x„,(7„) := (e^*^°(*-i(xi)...*'^"(a;„))e-^*^°). (40) 

Let us also introduce so-called truncated Wightman functions according 
to following relations: 

Wn,t{xi,ai,...,Xn,(Jn) = ^ JJ Vl^|p|,t(a;iP , (7^P , X^p ^ , (T^p ^ ) . (41) 

Here \ A\ is a number of elements of finite set A, P is a set of all decompositions 
TT = {Pi, i-|7r|} of the set {l,...,ri} into disjoint union of sets Pi, i-|7r|. 



14 



For each k = 1, |7r| we put by definition = {«f*, ...,i^p^^}, zf*" < < 
^ -Pk 

■■■<'\k\- 

Note that the property of enough fast asymptotic space factorization 
of correlations by definition means that truncated Wightman functions 
cxi, x„, o"„) tends to zero enough fast if their relative arguments 
= X2 — xi, ...^n-i = Xn — xi tcuds to infinity. 

Now let us introduce the Fourier transformations of Wightman functions 
as follows: 

Wn,t{pl,Cfl, ■.■,Pn,0-n) ■ 

= J d^x,,...,d^Xne'^^^P''''+---+''-P-^-^WnM,ai,.^^ (42) 

and 

^ltiPu(^i^ ...,p„,(T„) : 
= J d'x,,...,£xne'^''^P^'''+---+^-P-^-^Wl,{xi,a,,.^ (43) 
one can prove that 

W^n,t(Pl,(7i, ...,Pn,CTn) = ^ ]^ H^p|_t(p,P , (T^P , . . . , p^P ^ , a^P ^ ) , (44) 

TTGP Pe-TT 

where notations are such that as in f HTj) . 

The truncated Wightman functions W^^{xi, fJi, x„, cr^) are translation- 
invariant and the function of (enough) fast decay on relative arguments. 
Therefore their Fourier transforms W^^ could be represented as follows: 

Wn,tiPi,(^i, ...,Pn,o-„) = 5{pi(Ti + ... +p„cr„)>Vji(pi,ai, ...,p„,o-„), (45) 

where Wj^ is an enough smooth function of at most polynomial increment 
at infinity (with enough large number of derivatives). 
It is obvious that Vn = 1, 2, ... and Vt G M 

Wj(pi,ai,...,p„,a„) = e^('^^-(^^)+---+--(^"))*Vr„^o(Pi,^i,---,Pn,^n)- (46) 

Now, let us prove that for each n = 3,4,... W„,i(pi, cti, cr„) tends to 

zero as t — ^ — oo in the sense of distributions. This means that for each test 
function /(pi, ...,p„) G ^(M^") 

j d?pi...d^PnWl^{pi,ai,...,Pn,an)f{pu...,Pn) (47) 
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as t — i- — oo. It is obvious from remarks above that to prove this fact it is 
enough to prove that for each enough smooth function g{xi, x„) of enough 
fast decaying at infinity 

J (i3p,...c/3p^e*K-(fi)+-+-n-(P"))*<5(p^ai + ... +Pn(Tn)g{pi, .■.,Pn) ^ (48) 
as t — )■ — oo. 

Put by definition L = {(pi, ...,pn)\aipi + ...anPn = 0}. Note that for each 
n = 3,4, ... aiUi{p) + ... + anOJnip) is not identically equal to zero on L. This 
fact could pe proved very simply but we omit this proof. So the limit equality 
fHHj) could be directly obtained by using the stationary phase method. 

So, for each n = 3,4,... iy„ ^(pi, ai, cr„) — )■ as t — )■ — oo in the 

sense of distributions. From other hand, it is obvious that VFj^ = W^^q- 
result Vt e M 

Wn,t{xi,ai, ...,Xn,crn) = (49) 

if n is odd and 

lim Wn,t{xi,ai,...,Xn,(yn) = y^^Y\W2,o{xi^,ai^,Xj^,aj^), (50) 

reTTer 

if n is even (in the sense of distributions). Here T is a set of all decomposition 
T = {Ti, ...,T^} of the set {l,2,...,m} into disjoint union of pairs = 

But let us recall that 

lim Wn,tixi,ai,...,Xn,an) = {'iJ'''{xi)...Psi''"{xn))o. (51) 

t— >— oo 

in the sense of distributions. Therefore is finite and Gauss. Now let us 
formulate main theorem of present section 

Theorem 2. let (■) be some stationary translation invariant state on the 
CCR-algebra, commuting with the number particle operator. Suppose that 
(■) satisfy to the asymptotic property of (enough fast) space factorization of 
correlations. Then 

Pi-) = |tr((-)e-^) (52) 

for some T > and /i G M. Z is a statistical sum defined by equation (!23|) . 
In other words the state (■) is a Gibbsian state. 



16 



Proof. It follows from lemma 5 that 

(•) 

= (>S-^(0,-oo)(-)-S(0,-oo))o. (53) 

This imphes that 

(T^|(ti,xi)...^|(t„,a:„)) 
= (5-iT^±(ti,.xi),...,^^(t„,x„)5)o. (54) 

But one can prove (in the sense of formal series of perturbation theory) 
that the fact that (■) finite implies that (T\l'^(ti, a;i)...\E'^(t„, a;„)) is finite. 
Therefore the Keldysh diagram technique for {T^f^{ti, Xi)...^^(t„, Xn)) does 
not contain divergences. Therefore lemma 4 implies that 

{a+ik)a+ik'))o = {a{k)aik'))o = 0, 

{a{k)a^k'))o = S{k - k')^^ (55) 

e T — 1 

for some T > and e M. And it follows from lemma 2 that 

(•) = |tr((-)e-^). (56) 

Therefore theorem is proved. 

Remark. Note that we have supposed the following when we have proved 
the theorem 2 above. If the state (•) is a finite state (does not contain the 
divergences) then the divergences does not exist in each order in coupling 

constant A of the power decomposition of the state (■) according to pertur- 
bation theory. We have not been prove this statement but it is essential to 
suppose that this statement holds for all states of the physical interest. 

3 The Bardeen — Cooper — SchriefFer 
model on the language of generalized ther- 
modynamics. 

In this section we show how to formulate the Bardeen — Cooper — Schrieffer 
model of suoercondactivity on the language of the generalized thermodynam- 
ics. But at first let us recall several results on some model Hamiltonian in su- 
perconductive theory which belong to Bogoliubov — Zubarev — Tserkovikov 
[5,9]. 
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Let us consider a model dynamical system described by the following 
Hamiltonian: 

H^Yl - ^ E Kmf'Hatja.f,af,. (57) 

/ /',/ 

Here we use the following notations: 

f={p,s), -f={-p,-s). (58) 
p is a particle's momentum, s = ±1, 

2 

T{f) = 1^ 0. (59) 

A(/) = j£(/) if l^-H < A 

A(/) = Oif |^-//| > A 

e{{p,s))^s, (60) 

Qf, a J are usual fermion annihilation-creation operators. 
We will also consider the Hamiltonian H with some additional terms 
which are the sources of creation and annihilation of pairs. 

H.-H-'^Yl ^(/){«-/«/ + 
/ 

u>0. (61) 
Let C be an arbitrary complex number. Let us represent H as follows: 

H^Ho{C) + H^{C), (62) 

where 

Ho{C) = TifH^f -^Yl Kf){C*a.faf + Ca^atf} + ^-^V (63) 

/ 



and 



^i(^) = E - E Kf)a-faf - C) (64) 



/ / 
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Note that Hq{C) is a quadratic form on fermion annihilation-creation oper- 
ators. So it may be diagonalized by mens of canonical fermion Bogoliubov 
transformation. More precisely: 

f 



where 



v,^-%l- , ^^^^ (67) 

Note that u{—f) — u{f) and v{—f) — —v{f). Function u is a, real- valued 
function and function v is a complex- valued function, and + \v\^ — 1. 
These implies that the operators 

«/ — ttfUf + atfVf, 

— a'jUf + d-fV*^ (68) 

are the operators of the fermion type. 
Put by definition: 

J{C):^^Y.^{Ma_f-C\ (69) 

/ 

Now the expression for Hi{C) could be rewritten as follows: 

Hiic) = -Y^r{c)j{c). (70) 

Let r be an arbitrary (low bounded) Hamiltonian. Denote by (■)r,Tthe av- 
eraging with respect to the Gibbs distribution corresponding to the Hamil- 
tonian r and the temperature T: 

tr((-)e-^) 

v/r,T — ———rr- y'^) 

tr(e Tj 
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Let us consider the following equation: 

(J(C))ho(c),t = (72) 

on the variable C . This equation is called the equation of compensation. Note 
that Hq{C) is a quadratic form on creation-annihilation operators. Therefore 
the equation of compensation could be rewritten explicitly as follows [5]: 



2V Y v/A2(/)|C|2 + T2(/) 2T s y ) 

In the limit V ^ oo this equation will take the form: 



=th{ — }. (74) 



(27r)2y ^x\p)\C\^ + T\p) 2T 

It is well-known fact that last equation has a non-zero solution for all temper- 
ature T < Tcr for some critical temperatures T < T^- for some temperature 
Tcr- We will consider only the case when T < Tcr- Let Co be an arbitrary 
solution of the equation ( 1721) or equation ( 1731) . Note that if Co is a solution 
of equation (173|) then Coe*"^ is a solution of equation (!73l) too for an arbi- 
trary real number Lp. Let Cq is a positive solution of the equation ([Tj) or 
equation (!73|) . Now let us formulate main results of Bogoliubov, Zubarev, 
Tserkovnikov [9] concerned with the system under consideration. 

Let r be a Hamiltonian of some system contained in some volume V). 
Denote by FyIV^T) the free (Gibbsian) energy of this system corresponding 
to the temperature T. Then the following equalities hold 

Now, let us consider quasi-averages 

-< •••«/,•••«/.••• yH,T-= 

= .jj^^r^,}^"^ (■■■^l-^fs-)H.,T (76) 

Here {...a^,...af^...)H„,T are fermion creation-annihilation operators in Heisen- 
berg representation. The following equation holds [9,5]: 

^ ...a+....a/,... yH,T 
= lim (...a+ ...a/,...)^^(c^),T. (77) 
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Now let us consider the following Hamiltonian instead of Hy\ 

= -\T. A(/){a-/a/e^^ + a+at^e^^, (78) 
/ 

where is an arbitrary. We will have again the equality f l75|) . or more 
precisely: 



where we put now 

l^"^^n ^-^X-^fs-^Ht^T (80) 



and 

Ct := Coe^^. (81) 

It is necessary to note that the equation of compensation ( |73l) is equal to a 
condition of the minimum of free energy FiKc^^^^iy^ T) on C at fixed energy T) 
(see [5]). Recall that Fh^(^c) is a free energy of the dynamical system which 
is described by the Hamiltonian Hq{C) and contained in volume V at the 
temperature T. But the condition that Co corresponds to the minimum of 
free energy Fho{Co) is equivalent to the following condition: the entropy of our 
system Sho{c} achieve the maximum at the point C = Co (as a function of C) 
at fixed energy E. Last remark admit us to suppose that thermodynamics of 
our model could be described in terms of our generalized thermodynamics and 
the parameter C here plays the role of the observable values of commuting 
integrals of the system. 

Let us recall that one of the basics principles of Generalized thermody- 
namics states the following. The observable values of of the commuting first 
integrals (from microcanonical distribution) which may be realized in nature 
must corresponds to the maximum of the entropy (under fixed energy). 

Let us start to reformulate thermodynamics of our model in terms of 
the generalized thermodynamics. It follows from the basics principles of 
statistical mechanics that one can use the microcanonical distribution: 

Pe{C) := const6{H^{C) - E) (82) 
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instead of the canonical distribution 

H„(C) 

e T 

= , , go(C),- (S3) 

tr(e T ) 

the entropy corresponding to the distribution has the form: 

S{E, C) = tr5(//o(C) - E). (84) 

Note that in equation fl82l) and flH^ by symbol 6{x) we mean some regular- 
ization of (^-function (5-function A{x) defined as follow: 

A(x) = 0, if Ixl > — 

2 

A{x) = ^, if \x\ < ^ (85) 

where A is some small positive real number (asymptotical constant as y — )■ 
oo). 

Let us consider the following symplectic manifold ©2 with the canonically- 
conjugated variables Sr, Sj, ^r, $/ on it defined as follows. 

Variables Sr, Sj run the whole real line, but variables ^r, $/ are defined 
modulo the substitutions: 

$R I-)- $/ + 2?™, 
$/ I-)- $/ + 2?™, 

(86) 

where n G Z. The Poisson brackets of variables 5'^, Sj, ^r, $/ by definition 
has the form: 

{Sr,Sj) = {<^r,<I>i)=0, 
{SR,<!>j) = iSi,<^R)=0, 
{Sr,<^r) = {Sj,<!>j) = 1. 

(87) 

Let ©1 be a (non-commutative) phase space of the dynamical system, de- 
scribed by the Hamiltonian f l57|) . Let (3 be a direct product of the (non- 
commutative)spaces ©1 and &2, & = &i x ©2- Now let us consider the dy- 
namical system on the space © which described by Hamiltonian Hq{Sr + Sj). 
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It is evident that the dynamical variables Sn and Sj are the commuting first 
integrals with respect to the Hamiltonian Hq{Sr + iSj). The generalized 
micro-canonical distribution corresponding to the integral Ho{Sii + iSj) and 
the integrals Sr and Sj could be written as follows: 

S'^iE, C) := trg. j d<pR...dSiS{H^{SR + iSj) - E) 

X 6{Sr - ^C)6{Sj - 3C) = S{E, C) + c. (88) 

Here S{E,C) is defined by equation f l51|) . c is some real constant which 
does not depend on E and C and symbol trei means trace over the Hilbert 
space of our model dynamical system described by the Hamiltonian flFTI) . 
Recall that one of the basics principles of our generalized thermodynamics 
state that the observable value of C (which may be realized in nature) must 
corresponds to the maximum of S'^{E,C) (at fixed energy E). We have 
proved that this condition is equivalent to the equation of compensation 
f l72p . Therefore, we have formulated the thermodynamics of our model in 
terms of generalized thermodynamics. Note also that all nonzero solutions 
of equation of compensation form some one- dimensional manifold. So the 
entropy S'®(£',C) (at fixed energy) has one- dimensional plateau. 

4 On some new method of derivation of 
Boltzmann kinetic equation. 

In the present section we consider Bose gas with weak pair enough good 
interaction. For this system we derive usual kinetic Boltzmann equation in 
the limit of weak interaction by some new method. We consider a single mode 
of Bose gas (corresponding to some momentum k as an oscillator interacting 
with other modes as an oscillator interacting with thermostat. To study 
the evolution of this fixed mode interacting with other modes (considered as 
thermostat) we will use the so called statistical perturbation theory [10]. 

One of the interesting consequence from non-ergoding theorem states the 
following. To obtain an irreversible macroscopic evolution from the reversible 
microscopic evolution one need to take into account the behavior of the sys- 
tem at its boundary. In other words it is impossible to prove that the ther- 
malization takes place in the system if we neglect by the role which plays the 
environment of the system in this process. The derivation of kinetic equa- 
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tions presented in this section gives us a new point of view at the role which 
plays the environment of the system in the process of its thermalization. 

So let us consider a model dynamical system, described by the following 
Hamiltonian, 

H^Ho + \V, (89) 
where Hq is a free Hamiltonian 

Ho^ / (y ~ fi)a+{k)a{k)(fk (90) 

and 

Hi = ^ J d^Pid'^P2d^qid^q2v{pi,P2\qi,q2) 
xS{pi +P2-qi- q2)a'^{Pi)a'^{P2)a{q2)a{qi) (91) 

Here a+(A;), a{k) are usual Bose creation- annihilation operators (described 
in section 2), e R is a chemical potential A e R is a couphng constant. 
v{pi,p2\qi,q2) is a test function from the Schwartz space such that 

v{Pi,P2\qi,q2) =v*{qi,q2\pi,P2). (92) 

Note that the Hamiltonian which describes Bose gas with weak (enough 
good) is a Hamiltonian of such form. 

Now let us briefly describe the representation of quasi-particles in sta- 
tistical mechanics. We will use this representation below when we try to 
establish the form of the density operator for non-equilibrium state. 

Let us consider Bose Gas contained in macroscopic volume V. Let us 
formulate the main physical assumption which lays in the fundament of 
quasi-particles representation. The classiflcation of the energetic levels of 
interacting particle could be described in the same manner as a classification 
of energetic levels of the system of non interacting particles. Therefore the 
energetic levels of the system could be described by means of real valued 
pointwise positive function n{p) of momentum p G M'^. One can interpret 
the density function n{p) as a density of quasi-particles with momentum p. 
In other words n{;p)d?p is a number of particles with momenta from infinites- 
imally small volume (in momenta space) d^p with the center at p. Let E[rii\ 
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be an energy of an energetic level corresponding to the function n{p). Let 
e{p) be a functional derivation of E[n] with respect n{p), i.e. 



One can interpret s{p) as a dispersion low of quasi-particles. We have sup- 
posed that the classification of energetic levels of interacting system is the 
same as a classification of energetic levels of non-interacting system. There- 
fore the entropy of the system corresponding to the density function n{p) has 
the form: 

S[n\ = J (fp[{l + n{p)) ln(l + n{p)) - n{p) Inn(p)]. (94) 

Note that the dynamics of (isolated) interacting system differs from the dy- 
namics of non-interacting systems only in the following/ The dispersion low 
for e{p, [ri\) for interacting system depends of the density function n. 
Therefore for all piece- wise continuous functions /(p) the quantities 

K[f] = J h{p)f{p)d^p (95) 

are commuting integrals of motion. Here, by definition 

h{p) := a~^{p)a{p) (96) 

and a'^{p) and a{p) are formally defined operators of quasi-particles, satis- 
fying to the Bose-Einstein statistics. 

Note that last observation is the main idea of the proof of non-ergodic 
theorem. More precisely for the systems which contains of finite number 
of particles (quantum fields) wave operator W (defined in section 2) maps 
dynamics of free particles into dynamics of interacting particles. Therefore 
it is essential try to find the construction analogues to the construction of 
wave operators to build the quasi-particles representation. But the series of 
perturbation theory obtained by such a way contain (secular) divergences. 
One of the authors comes to non-ergodic theorem when tried to renormalize 
such divergences. 

Let {Oi\i = 1,2, 3...} is a decomposition of M.^ into a set of domains such 
that Vi ^ j OiH Oj C dOi U dOj , where dU means the boundary of the set 
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U. Suppose that yi — 1,2,3...N and Vp,p' e the difference p — p' is very 
small. Vi = 1, 2, 3... put by definition 

Ki^ J n{p)d^p (97) 

Oi 

Now let us find the form of distribution function n{p) corresponding to some 
(generally) equilibrium state (in the sense of generalized micro canonical dis- 
tribution) which corresponds to fixed values K'^ of integrals Ki, i — 1,2,3, .... 

It is evident that (generally) equilibrium distribution function n(p) cor- 
responds to (relative) maximum of the entropy S{n) under the following 
constrains: 

K,[n] := J d?pn{p) = K[, 

1 = 1,2,3.... (98) 



and 



E[n] = E. (99) 



According to the method of undefined lagrange multipliers we reduce this 
problem to the solution of the following system of equations 

^{-EG.l(l+e-)>n(l + ^)-g-lng-] 



i=l 



N 



J2f^iNi + E[n]} = 0, (100) 



i=l 



where 



:= j d'p (101) 



and 



= y n{p)£p, 

Oi 

i = 1,2,3... (102) 
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This system of equation could be immediately solved and one finds 

^iP) = (103) 

e T — 1 

where 

Ai(p) = e M, Hp e Oi (104) 

and T some positive number 

In other words if one moves from generalized microcanonical distribu- 
tion to the generalized microcanonical distribution this is equivalent to the 
following. The chemical potential /i becomes to the function of momenta p. 

Note that in spite of volume V is macroscopic it is finite. Therefore the 
set of all possible momenta is discrete {pi\i — 1,2,3...}. Let by definition Hj 
be a number of particles with momentum pi. 

Wi = 1,2,3, ... the cells Oi corresponds to the approximately equal fre- 
quencies. Therefore each two states corresponding to the same numbers 
Ni, i — 1,2,3,.... corresponds to approximately the same energy. There- 
fore all the states corresponding to the same numbers i — 1, 2, 3, ... have 
approximately the same probabilities. 

One can easily find the form of the distribution function pi{ni) by using 
the method of the most probable distribution. Omitting the well-known 
reasoning one finds 

Piirii) = f , (105) 

where 7,^ is a real positive number such that 

oo 

Y.Pi{n)^\. (106) 

n=0 

Therefore the (generalized) equilibrium state of quasiparticles could be 
described by means the Gauss state on the algebra of canonical commutative 
relations generated by Q;(p), ol^{p) of the form 

(a(;>)a(p')) = («+(p)a+(p')) = 0, 
(a+(p)a(p')) = <^(/>-p')(l + ^(/>)), 

n{p) = S{p-p')^^;^^ , (107) 

e T — 1 
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fxip) = fXi, Hp e Oi. (108) 

So, let us consider a single mode of momenta ki interacting with other 
modes of Bose gas as with a thermostat. To analyze this system we will use 
so-called statistical perturbation theory [10]. Bogoliubov [10] has considered 
the following system. He has considered the Hamiltonian system of n degrees 
of freedom described by the Hamiltonian H{p, q) in the field of external force 
ef{t), £ e R. So, total Hamiltonian of the system has the form 

T = H + ef{t)P, (109) 

where P{p, q) is some function of canonically-conjugated of momenta and 
coordinates. We also suppose that external force f{t) is the function of time 
t of the form: 

/(t) = ^a^cos(z^t + ^^). (110) 

Here the phases t/?^ are independent random quantities uniformly distributed 
on circle. We also suppose that the spectrum of frequencies is almost con- 
tinuous so we can the sums of the form 

(111) 

V 

for continuous functions F{y) replace by the integral 

o 

F{v)l[v)dv. (112) 



Let Dt denote the probability density of momenta and coordinates at 
fixed time moment t under the condition that phases t^)^ take some defined 
values. The probability density of momenta and coordinates in usual sense 
(for arbitrary values of phases ip^) could be obtained from Dt by means of 
averaging over the phases ^fy. In other words 

D^ = po. (113) 

The evolution of Dt in times could be described by means well known Liou- 
ville equation 

dD 

-^ = {H,Dt) + efit){P,Dt), (114) 
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where {A, B) denotes usual Poisson Brackets defined by the formula 



(A,B).±iff-ffy ,115) 

dqi dpi dpi dqi 

Let us introduce the following one-parameter group of operators Tf acting on 
dynamical variables according to the formula 

TtF{p,q) = F{pt,qt), (116) 

where [pt, qt) are the solution of canonical Hamiltonian equations 

dpt _ dH{pt,qt) 
dt dqt ' 

dqt ^ dH{pt,qt) 
dt dpt 

corresponding to the following initial conditions 



(117) 



Po = P, 

Po = P- (118) 

In these assumptions and denotations in the limit of small coupling constant 
e it is derived [Bogoliubov, 1945] the following equation on pt 



t 

dpt 
dt 



{H,pt) + e^ J A{t-T){P,{Tt^r{P,Pr)))dT, (119) 



where 

+ 00 



A(r) = ^ J cos{i^T)dT. (120) 



The equation flll9p could be essentially generalized to the case of quantum 
mechanics. It is enough to replace the Poisson brackets by the commutator. 

Note that the random force f{t) here is an example of random processes. 
Recall that the random process is the set (T,ft), where T = (5, S,P) is a 
Kolmogorov triple {S is some set, S is a a- algebra, P is a probability measure 
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on S defined on sets from S). ft is a map wliicli to an arbitrary f G M assigns 
a random quantity i.e. a measurable (witli respect to S) function on S. Note 
that the equation (11191) could be rewritten as follows 



dpi 

at 



[H, pt) j JWi^iP, m^AP, Pt)))dT, (121) 







where ^ means the mathematical expectation of ^ i.e. 

(122) 

s 

Note that how it follows from the derivation of equation (I119P [Bogoliubov, 
194] it is not necessary that f{t) has the form described previously. It is 
enough that f{t) has the form 

where fu are random quantities such that Wp, ly, 

^ = ij=0, (124) 



Ufu = f;f* = 0, (125) 

and 



= (126) 

The spectrum of frequencies {v} is supposed almost continuous such that 
one can replace the sums 

^J,F(z.) (127) 

V 

for continuous functions F^v) by the integral 

+ 00 

j I{v)F{v)dv. (128) 
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Now note that the generahzation of a Kolmogorov triple to noncommutative 
(quantum) case is a pair consisting of an unital ^-algebra A and a state (■) 
on it. By definition state is a normahzed positively defined linear functional 
on the algebra. Normalization condition means that (1) = 1, where 1 is a 
unite of A. Positively definiteness means that "^a E A (aa*) > 0. 

The algebra A plays the role of the space of measurable functions on 
(nonexisting) noncommutative space. The state (■) plays the role of inte- 
gration on this state. The involution * plays the role of pointwise complex 
conjugation. 

By definition the noncommutative random process is a map which to an 
arbitrary t G M assigns an element at of A. We say that noncommutative 
random process {at\t G M} is a real-valued process if Vt G M = a*. 

The equation flllQp could be generalized to the case of noncommutative 
(quantum) systems. More precisely let total Hamiltonian of the system has 
the form 

H = Ho + s{fit)P + f+it)P+), (129) 
where f{t) is a noncommutative random process of the form 

(/.) = (/;) = 0, (131) 

and V/i, V 

(fju) = {f;ft) = 0, (132) 



iUt) = Ju6,u. (133) 

The spectrum of external force f(t) is assumed almost continuous such that 
we can replace the sums 

V 

Y,F{v)J, (134) 
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by the integrals 

+00 

F{u)I{ij)du, 





+00 



F{v)J{v)dv (135) 







respectively. In this case, instead of f lllQp we will have 

t 

+e' lmt)P\f{r)T^_,{[P,Pt])]) 




t 



+e' l{[fmr{r)T^-t{[P\pt])]). (136) 


We suppose that our system is contained in some cube of volume V with 
edges parallel to the coordinate axis. We suppose that the length od the 
edges of this cube is equal to L. Therefore V = L^. Put by definition 

^={*=(t..,.y|^,^,J^.Z}. (137) 

We need to use the discrete momentum representation for finite volume in- 
stead of continuous momentum representation for infinite volume. Instead 
of the Hamiltonian fl89|90t WT\i we need to use the following Hamiltonian in 
finite volume: 



where 



and 



= + hX, (138) 
= ^uj{k)a^{k)a{k), (139) 



u{k) = y. (140) 
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Hi = 2^ ^{pi+P2- qi- q2)v{pi,P2\qi,q2) 



Pl:P2;«l,«2€J? 



xa+a+Og^agi, (141) 



where 



A{p) = 1, if p = 

A(p) = 0, if p 7^ 0, (142) 

and and are creation and annihilation operators in discrete momentum 

representation. 

Now let us consider a single mode corresponding to some fixed momentum 
k E A. Let us study how the state of this mode changes while time changes 
from t to t + dt for infinitely small dt. To study how the density matrix 
Pt of this mode changes at time interval [t, t + dt] one can suppose that 
the thermostat evolves at this interval according to free dynamics. So one 
can suppose that the mode corresponding to the momentum k affected by 
external (noncommutative) force defined by this thermostat. In other words 
we must take 

r^H + X{f{t)P + r{t)P+), (143) 

where, by definition 

H ^u(k)a+(k)a(k), 

P = afe. (144) 

External forces f{t) are the elements of the algebra A generated by operators 
Op , Op, p ^ k. External force f{t) has the form 

(27r)'^ 

f{t) = —^ ^ A(pi +P2 - ?i - A;)v(pi,p2|gi,g2) 



Pl,P2,?l6il\{fc} 



To obtain the differential equation for the density matrix for our fixed mode 
we will average over the state (■) on the algebra A which is a Gauss state 
and is defined by the following its pair correlator 



a. 



pOp') = (0+0+) = 0, 



(a+ap,) = A(p-p')n(p). (146) 
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Note that external force f{t) has the form 

V 

where 

/- = ^ E A(pi+p2-gi-fc) 

pi.P2,giG^\{fc} 

xA(a;(pi) +a;(p2) -^(gi) - i^)^^(Pi,P2|gi, 92)0^1 Op^Sa- (148) 

Note that the averaged pair products of j^^ f* satisfy to the equations 
fll32lll33p . Let pt be a density matrix corresponding to our fixed single mode 
with momentum k. Put by definition (■) = tr(-pt). 

Suppose that (■) be a Gauss state on the algebra A. Equation (11361) 
implies that Vt > the state {■)t+dt is a Gauss state if the state {■)t is a 
Gauss state for infinitely small positive dt. Therefore if {■)t is a gauss state 
for all positive t if is a Gauss state (in the limit of weak pair interaction). 
Let be an operator of the particle number corresponding to the 

momentum p. Put by definition nt{p) = {np)t- Equation (I136p implies that 
(in the limit of weak pair interaction) 

t 

jntik) = -X'l{{[atrit),aJir)]),)e-^'^'~'-^dr 



t 



= -A^ f{ntmnt)f{r)) - {l + n,{k))U{r)nt))}e'-^^^^'-^^dT 
Jo 

+X' A(l + n,(fc))(/(t)r(r)) -n,(A;)(r(r)/(t))}e--('=)(*--)f/r 
Jo 

t 

= -XWit) j (r(0)/(r))e-^-«-(ir 
-t 

t 

+\\l + n,{t)) J{f{0)r{T))e'-'^'^^dT. (149) 
-t 
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Therefore 

d 



= 775— 22 \'"iPuP2\<li,k)\nt{qi)nt{k){l + nt{pi)){l + nt{p2)) 



y2 

pi,P2,qieJ?\{fe} 



X 



A(pi +P2-qi-k) J e*(-(Pi)+-(P2)— (</i)— fe))r^^ 
-t 

A^(27r)^ 

..2 XI b(Pi>P2|gi,A;)|(l + nt(5i))(l + nt(A;))nt(pi)nt(/)2) 



pi,P2,«ieJ?\{fe} 



t 



-t 

But volume V is macroscopic. So we can replace the sums over ^ by integrals. 
In result 

jnt{k) = (27r)A^ / d'^pid^p2d\iS{pi + p2 - Qi - ^2) 

5(a;(pi) + a;(p2) - a;(gi) - x {(1 + nt{qi)){l + nt{k))nt{pi)nt{p2) 

-nt{qi)nt{k){l + nt{pi)){l + nt{p2))} (151) 

So we have obtain the usual kinetic equation for Bose gas in the limit of weak 
interaction. 

5 Integrals for non-linear partial differential 
equations. 

Soliteiry perturbations from generalized thermodynamics. Let us 

consider some system of statistical mechanics such that there second-order 
phase transition takes place at some temperature Tcr- Suppose that this 
transition concerning with some symmetry breaking. As an example of such 
system one can take some superfluid or superconductive system. Let (•) be 
(generally) equilibrium translation-invariant stationary state corresponding 
to the temperature T < T„.. Suppose that (P) = 0, where (P) is a mo- 
mentum operator. Let /C be some inertial reference system which moves 
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with velocity v. Let U C be some domains of microscopical sizes tightly 
connected with reference system /C. Domain U by definition is macroscopic. 
Therefore, according yo the non-ergodic theorem [1] there exist integrals 
Ki, Kn in involution for the part of our system contained in U. Accord- 
ing to our generalized thermodynamics, the distribution function for the part 
of our system, contained in U has the form: 

N 

p{x) = const5{H{x) -E)W 5{K'^{x) - K'^). (152) 

1=1 

Now let us show, that the entropy of the system contained in U : 

N N 

S{E, K[, O = In / drj{H{x) -E)l[ 5{H{x) - E)]^ 5{Ki{x) - K'i)15S) 

i=l i=l 

achieve a maximum at some D > dimensional manifold M as a function 
of integrals Ki, ...,Kn at fixed energy E. Indeed suppose that the entropy 
achieve a maximum (at fixed energy) at isolated point. In this case the (gen- 
eralized) microcanonical distribution gives us the same results as standard 
microcanonical distribution 

p{x) = constS{H{x) - E). (154) 

But for this distribution all the symmetries of the Hamiltonian H{x) are 
preserved. Therefore the entropy S{E, K[, K'j^) achieve a maximum at 
some M > dimensional manifold as function of K[, K'j^. 

Therefore we can chose the observable values K[, of the integrals 

Ki,...,Kn such that (the intensive) properties of the system will be quali- 
tatively different from the analogues properties of other part of the system 
(outside of U) and values K[,...,K'j^ corresponds to the maximum of the 
entropy. According to our generalized thermodynamics this state is stable. 
For example, for superfiuid gases, the the phase of condensate wave function 
inside U may differ from the phase of condensate wave function outside of 
U. If we return into initial (static) reference system we obtain a non-trivial 
perturbation of the state (•) which moves with constant velocity without 
dissipation. Such perturbation we will call thermodynamical solitons. 

An example of thermodynamical soliton. The Abrikosov Vor- 
tices. As an example of thermodynamical soliton we consider the situation 
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which arise in the superfluidity theory if one studies the condensate wave 
function [11]. 

Let \E'(x) and be secondary quantized wave function and its 

hermitian-conjugated. Let S and S"'"(a;) be the parts of the function ^'(a;) 
and ^"'"(a;) which change by one number of condensate particles and preserve 
others quantum numbers. In other words, by deflnition: 

S(a;)|A^ + l,m) = E(x)|iV,m), 
E+{x)\N,m) =J:*{x)\N +l,m) (155) 

Here 

E{x) -.^ {N,m\E{x)\N+l,m), (156) 

and E*{x) is a complex-conjugated to E{x) function. TV is a number of 
condensate particle m is a number of other quantum numbers of the system. 
Let us represent ^(x) and ^(x) as follows: 

^{x)^E{x) + ^'{x), 
*+(x) = E+{x) + ^'+{x). (157) 

In thermodynamical limit N ^ oo the difference between the states | A^, m) 
and \N + l,m) disappears at all. Therefore the operators E{x) and S+(x) 
becomes the operators which commutes to each other and commutes with 
\i''(x) and \l/'+(a;). In other words E{x) and E'^{x) becomes the classical 
variables. Now let us write the (partial) differential equation on condensate 
wave function E{x). Let us consider a weakly non- ideal Bose gas at absolute 
zero of temperature. Almost all particles in such a gas are in condensate state. 
In terms of secondary-quantized wave function this means that the over- 
condensate part \l/'(a;) is smafl with respect to the condensate wave function 
E{x). If we neglect by small over-condensate part ^'(x) the condensate wave 
function will satisfy to the same Schrodinger equation as the equation which 
take place for the operator ^(x). I.e. 

+E{x,t) J \E{x',t'\^U{x-x')(fx'. (158) 

Here U{x) is an interaction potential /i — nUo, n is a density of particles in 
gas and Uq J U{x)d^x. This equation is called the Pitaevskii equation. 
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We will show now that the Pitaevskii equation will precise in the limit when 
the interaction constant tends to zero. Let us suppose that the interaction 
potential U {x) depends on small positive parameter A as follows 



U{x\X) := XX^/^U{xv^). 



(159) 



For example if U{x) = S{x) then U{x,X) = XS{x). Let us introduce new 
condensate wave function: 



In the limit A — )■ 0, S'(a;,t) will satisfy precisely to the Pitaevskii equation 
(11581) . The fact that the Pitaevskii equation will be satisfied precisely will 
follows from the fact that all particles will be in condensate in limit A — )■ 
and the over-condensate wave function will be equal to zero in the limit 
A ^ 0. 

The (partial) differential equation for the order parameter of 
the system has a lot of non trivial first integrals. In this subsection 
we consider only super fluid systems. These systems are described by the 
condensate wave function 'E{x) defined on M^. The condensate wave function 
play the role of order parameter for the system. In the limit A — )■ the 
condensate wave function S(x) and its complex conjugated H*(x) are the 
canonically-conjugated coordinates on some (infinitely-dimensional) phase 
space. The Poisson brackets of H(x) and 'E*{x) has the form: 



The evolution of the condensate wave function ^(a:) is a Hamiltonian evo- 
lution with respect to the Poisson brackets (I16ip . The corresponding 
Hamiltonian (for an arbitrary temperature) is F(S(x), S"^(a;))|T), where 
F(E{x),E~^{x))\T) is a free (Gibbsian) energy of the system under the con- 
dition of fixed condensate wave function S(x). Now we will introduce the 
assumption of so called asymptotically completeness and prove that the equa- 
tion for order parameter has a lot of first integrals in involution under this 
assumption. 

Asymptotical completeness. Let S(x, t) be a solution of the Pitaevskii 
equation (I158p which is localized (in some essential non rigorous sense) in 
some finite domain of the space at t = 0. We ssaume that this solution splits 
as t — 7- oo into the set of separated thermodynamical solitons described in 
the beginning of this section. 



s'(x,t) = j:{xVx,tx). 



(160) 



(H(x),H*(x')) = 6{x-x'). 



(161) 
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So let us suppose that the asymptotical completeness takes place. 
In this case first integrals for the system, defined by the Hamiltonian 
5+(2;))|r) (for the temperature T = 0) could be described as fol- 
lows. 

Let L be a cubic lattice in R^. Let {Ci\ i — 1,2,3, ...} be the set of all 
elementary cubes of L. We suppose that the sizes of cubes Q are macroscopic. 
We will also denote by Q (Vi = 1,2,3,...) the part of our system contained in 
the cube Cj. Vi = 1,2, ... system Ci is macroscopic. So Let {Kf}, j = N he 
some set of its first integrals in involution. We have shown [2] that number 
N could be chosen as large as needed if each system Q is enough large. Let 
Va e r„ : r3 jg ^ ^^^^^p defined as follows: 

Ta-.x^x + a. (162) 

Va e put by definition C(a) = T^C. 

Let a be an arbitrary vector from M^. Denote by (a), j = 1, 2, 3, ... the 
integrals of the system Q(a) obtained from Kf by translation on vector a in 
obvious sense. 

Let Va e ^H^) — {Kf{a))s,E*,T=o- Here {■)e,e*,t by definition is an 
averaging by relative Gibbs distribution corresponding to fixed S(x) and S* 
and the temperature T. 

Now let S(.x) be some localized in the space function. Let S(a:,t) be 
a solution of Pitaevskii equation such that S(x, 0) = S(x). Recall that 
Pitaevskii equations are Hamiltonian equations with respect to Hamiltonian 
F(E,(x), S"*" (x) ) I r) . According to the asymptotic completeness assumption if 
chose enough large time t we can assume that S(.t, t) is a set of space sepa- 
rated thermodynamical solitons such that for each i = 1,2,3... each cube Cj 
has nontrivial intersection at most with one of such solitons. Vj = 1,2, ...A'^ 
put by definition: 

/oo 
d'aJ2K^M^M,E*(x,t)). (163) 

It is evident that this definition do not depend on t (if t is enough large) and 
/Cj(S, S*) are integrals of motion. 

Theorem. Integrals /Cj(S, S*), j — 1, 2, N are in involution i.e. com- 
mutes to each other. 
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Proof. Suppose that the time t is chosen as above. We have: 



(/c„(h,e:*),/c„(h,s*)) 



J JCixCi i,j=l,2,... 

J JcixCi J ilR3xR3 



X 



X 



+ / / (fdidrd2 / / d'^xid'^X2 

J icixCi J Jr3xR3 

^ ' SE{X2) ("0,^(^2)). (164) 



Therefore to prove that the Poisson Bracket (/Cm(S, S*), /C„(S, S*)) is equal 
to zero it is enough to prove that the foUowing equahties holds: 

z:rr' = ,r 7' r ' = (i65) 

5E{x) 6E*{x) ^ ' 

for all i = 1,2,3..., j ~ 1,...,N, a,x e M^. Let us prove second of these 
equalities. Second of these equalities could be proved by analogy. 

Fix a cube Q(a). Suppose that the soliton (at most one) which has 
nontrivial intersection with Cj(a) has zero velocity. In the opposite case we 
can chose an inertial system such that the soliton has zero velocity in this 
new system. Then we just perform all reasoning below in this moving system. 

The (non-commutative) phase space of the system Ci{d) could be repre- 
sented as a direct product of two non-commutative phase spaces ©i x 62 
where the phase space ©i is a classical space and E{x) and E*{x) are the 
set of canonically-conjugated coordinates on it. And the space 62 roughly 
speaking corresponds to over-condansate particles. Let be a Hilbert space 
corresponding &i. Hamiltonian H and the integrals Kj of the system Ci{d) 
are the functionals of E{x) and E*{x) which take values in the space of op- 
erators acting in H. We have 

,^ _ tr(i^^(S,S-)exp(-^)) 

tr(exp(-^)) ■ ^ ^ 
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Here trace is taken over the space 'Hi{a). We have 

^^L(S,S) _ tr{^^ffpexp(-^)} 
<^S(x) tr(cxp(-M^)) 

{trX^(S,S*)exp(-. ^y))q|f)exp(.^y))exp(-^)} 



tr(exp(-^(^)) 



tr{K„(S, -)exp(-^)} tr{Mexp(-^)} 



■167) 



tr(exp(-^^)) tr(«^exp(-^) 

Let us show that each of three term in right hand side of this equation is 
equal to zero. Let us consider the third term. It is evident that this term 
is proportional to S(x,t). But the soliton which has nontrivial intersection 
with Ci{a) moves with zero velocity. So, this term is equal to zero. Let us 
consider first term: 

tr{ '^-jf;f ) exp(-^^^%g^)} 

Ijl^J . (168) 

tr(exp(-^(^)) 

This term is proportional to (H'(a;,t)), where dash means the derivative cor- 
responding to the flow, which generates by ir„i(S, S*). The averaging here 
is an averaging over generalized microcanonical distribution. The flow gen- 
erated by Xm(S,5*) preserve this microcanonical distribution. Therefore 

(S'(x,t)) =0. 

At last let us consider second term. Let us represent it as follows: 
^1 K ,^ JtrX^(S,S-)[g(S,S-),...[g(S,S-),q|f)]...]exp(-^)} ^^^^ 



n=0 



Note that the iterated commutator here remains flnite if the volume of Ci{a) 
tends to infinity. Note that the averaging in last formula is an averaging over 
the generalized microcanonical distribution: 



N 



p = const(^(if - ^) n ^(^i - ^j) (1 ''0) 
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In result in each term of series fll69p the integral Rmfn^E, S*) could be re- 
placed by its observable value K'^. In result the second term is equal to 



One can prove that this term is equal to zero by the same method as has 
been used for third term. Therefore the theorem is proved. 

Programm — Hypothesis. So we have proved that the equation for 
condensate wave function (hamiltonian with respect to F(H,S*|T)) admit a 
lot of independent first integrals in involution (in limit A = 0, T = 0). In the 
limit A = 0, T = the condensate wave function (after suitable rescaling of 
space coordinates and time) satisfies to Pitaevskii equation. Pitaevskii equa- 
tion is a non linear Schrodinger equation. In ID case the nonlinear Pitaevskii 
equation is an completely integrable system. But we have proved that in 3Z) 
case nonlinear Schrodinger equation (under the assumption of asymptotical 
completeness) admit a lot of independent integrals in involution. Therefore 
there arise the following program-hypothesis of constructing (not necessary 
completely) integrable systems. One can take some system of statistical me- 
chanics such that there exists a phase transition of second order at some 
temperature concerning with some symmetry breaking. Then one may find 
some small parameter rj for this system and write some (partial) differential 
equation for the order parameter in the limit A = by some exact asymptot- 
ical method. Obtained equation will be a hamiltonian system which admit 
a lot of independent commuting first integrals in involution. 

6 Conclusion. 

In the present paper we have considered some physical examples to illustrate 
basics principles of the generalized thermodynamics, developed in [2]. I am 
very grateful to A.V. Koshelkin for valuable critical comments on this article 
and very useful discussions. 
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